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ABSTRACT 

In recent years, a lot of work was done that has revealed some very interesting prop¬ 
erties of neutron stars. One can relate the first few multipole moments of a neutron 
star, or quantities that can be derived from them, with relations that are indepen¬ 
dent of the equation of state (EoS). This is a very significant result that has great 
implications for the description of neutron stars and in particular for the description 
of the spacetime around them. Additionally, it was recently shown that there is a four 
parameter analytic spacetime, known as the two-soliton spacetime, which can accu¬ 
rately capture the properties of the geometry around neutron stars. This allows for 
the possibility of describing in a unified formalism the astrophysically relevant prop¬ 
erties of the spacetime around a neutron star independently of the particulars of the 
EoS for the matter of the star. More precisely, the description of these astrophysical 
properties is done using an EoS omniscient spacetime that can describe the exterior 
of any neutron star. In the present work we investigate properties such as the loca¬ 
tion of the innermost stable circular orbit Risco (or the surface of the star when the 
latter overcomes the former), the various frequencies of perturbed circular equatorial 
geodesics, the efficiency of an accretion disc, its temperature distribution, and other 
properties associated with the emitted radiation from the disc, in a way that holds for 
all possible choices of a realistic EoS for the neutron star. Furthermore, we provide 
proof of principle that if one were to measure the right combinations of pairs of these 
properties, with the additional knowledge of the mass of the neutron star, one could 
determine the EoS of the star. 

Key words: gravitation - stars: neutron - equation of state - X-rays:binaries - 
accretion discs - methods: analytical. 


1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) and the neutron stars 
(NSs) that reside in them are sources of some of the most 
interesting astrophysical phenomena. They also provide the 
opportunity to study the most extreme physics in nature. 
The processes involved in these astrophysical objects as far 
as the NSs are concerned, are related both to strong gravity 
effects and effects associated with the properties of matter 
in densities that are above the nuclear density. Therefore 
these systems could be used as tools for broadening both 
our understanding of gravity, by testing the predictions of 
the established theory of General Relativity (GR), and our 
understanding of the equation of state (EoS) of matter for 
densities as high as the ones that can be found at the centres 
of NSs. 

Thus far, the uncertainty in the properties of the EoS on 
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the one hand and the fact that there are a lot of physical pro¬ 
cesses that are taking place in the interior and around NSs, 
has made them “dirty” systems for extracting clean signals 
for gravity. This picture partially changed recently with the 
discovery of the universal behaviour that NSs exhibit with 
respect to their structure, i.e., behaviour that is indepen¬ 
dent of the EoS used for the nuclear matter inside the NS. 
It was shown initially by Yagi & Yunes (2013b,a) that for 
slowly rotating NSs there exist EoS-independent relations 
between the normalized moment of inertia I = I /M 3 (where 
M is the mass), the normalized or reduced quadrupole 
Q = —QKM^j 2 ) 1 (where j = J/M 2 is the spin parame¬ 
ter and J is the angular momentum), and the normalized 
Love number A = A/M 5 . Pappas & Apostolatos (2014) gen¬ 
eralised the universality of the I — Q relations for arbitrar- 


1 This quantity will also be referred to as quadrupolar deforma- 
bility. 
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ily rotating NSs 2 and further showed that a more funda¬ 
mental universality holds, related to the spacetime proper¬ 
ties, that connects the first four multipole moments in an 
EoS-independent way. In their work, Pappas & Apostolatos 
(2014) additionally showed that if one were to measure the 
first three multipole moments of the spacetime around a 
NS, then one could distinguish an EoS for the star from the 
known realistic EoSs. The universal relations for the multi¬ 
pole moments in GR have been extended up to the mass hex- 
adecapole, Mu, which is the fifth non-zero moment, and have 
been shown to describe apart from NSs, quark stars as well 
(see the work by Yagi et al. (2014) and also see the work by 
Stein et al. (2014) for a generalised Newtonian treatment). 

For the LMXBs, there are several properties (or observ¬ 
ables in some cases) that can be associated to the properties 
of the spacetime and more specifically to the properties of 
geodesics around the compact object that is part of such a 
system. 

Quasi-periodic oscillations (QPOs) of the X-ray flux ob¬ 
served from LMXBs is one example (for a review see Lamb 
(2003); van der Klis (2006)), where one of the proposed 
mechanisms for producing this behaviour is the relativistic 
precession of geodesics (see Stella (2001)). This mechanism 
for explaining QPOs is not the only one that has been pro¬ 
posed, but it is the most straightforwardly related to prop¬ 
erties of geodesic motion. In the context of the relativistic 
precession model, one assumes that fluid elements of an ac¬ 
cretion disc follow geodesic or almost geodesic orbits which 
are almost circular and have some inclination with respect 
to the equatorial plane (an assumption that can be consid¬ 
ered valid for radiatively efficient accretion discs). General 
relativity predicts that these orbits will generally display a 
precession of the periastron of the orbit and a precession of 
the orbital plane, or of the nodes of that orbit. Therefore, 
there are several frequencies associated to these orbits, the 
orbital frequency of the circular motion Q, the radial oscil¬ 
lation of the fluid element along its orbit that has frequency 
k p and the corresponding precession frequency = Q — k p 
(periastron precession), and finally the vertical oscillation of 
the fluid element with respect to the equatorial plane which 
has frequency k z and the corresponding precession frequency 
Q z = Q — k z (nodal precession). Here, the term “frequency” 
refers to the angular velocity that is equal to 2ni/, where v 
is the frequency expressed as cycles in the unit of time. 

Another property that can be associated to the geom¬ 
etry of the spacetime around a compact object that is sur¬ 
rounded by an accretion disc, is the total luminosity radiated 
from the disc. This is a property that is associated to the 
efficiency r/ = 1 — Eisco of the disc, where -Eisco is the en¬ 
ergy per unit mass of a particle that has a circular orbit at 
the location of the innermost stable circular orbit (ISCO), 
which is essentially the total energy that a particle of unit 
mass loses as it moves from infinity up to the ISCO. If one 
assumes a thin radiatively efficient accretion disc, then a 
fluid element of unit mass in the disc will have radiated by 
the time it plunges from the ISCO to the compact object, a 
total amount of energy that is equal to 77 times its mass. 

Along the same lines, the temperature distribution of a 


2 Almost simultaneously Chakrabarti et al. (2014) came out with 
the same results and some further extensions. 


thin radiatively efficient disc can be associated to the varia¬ 
tion of E of circular geodesics with the radial distance from 
the central object. Therefore, the emmited spectrum of the 
disc can be associated to the geometry of the spacetime. 
This type of properties of an accretion disc and their connec¬ 
tion to the background geometry have been used in the past 
to determine the rotation of black holes (see for example 
Shafee et al. (2006) and references therein), where a tech¬ 
nique called, continuum-fitting method, is applied in order 
to determine the spin of the black hole from the luminosity 
of the disc. 

1.1 Executive summary 

This work is relatively lengthy and there are several different 
ideas that are being discussed. Therefore, in order to help 
the reader, we provide here a brief summary of the topics 
discussed. 

In what follows we give a universal, in the sense of EoS- 
independent or more precisely not EoS-specific, description 
of the spacetime around NSs and consequently a univer¬ 
sal and EoS-independent description of various astrophysi- 
cal phenomena that take place around these NSs. To have 
this EoS-independent description of the spacetime we use 
the analytic two-soliton spacetime, that has been shown to 
be an accurate representation of the spacetime around NSs 
(see work by Pappas & Apostolatos (2013)), in conjunction 
with the universal three-hair relations found for NSs in GR. 
This results in a description of the spacetime in terms of 
the mass and two more parameters, the spin parameter j 
and the quadrupolar deformability parameter a = Q. The 
mass essentially introduces only a scale and can therefore 
be removed from the various observables by applying the 
appropriate reduction. Therefore we end up with a para¬ 
metric description of the various mass-reduced observables 
in terms of only two parameters, j and a. The universal 
spacetime description is supplemented by a discussion of two 
more approximately EoS-independent relations that give the 
mass-reduced equatorial radius and the mass-reduced rota¬ 
tion frequency of NSs, in terms of the parameters j and a. 

The aforementioned observables have been separated 
into two categories, the observables that are related to the 
characteristics of geodesic motion on and around the equa¬ 
torial plane, which are mainly the various frequencies of the 
orbital motion, and the observables that are related to the 
radiation emitted by an accretion disc that forms around a 
NS that accretes matter. 

The observables discussed here, that belong to the first 
category, are the orbital and the nodal precession frequency 
of circular orbits at the location of the innermost available 
orbit around the NS. That orbit can either be the ISCO 
or the circular orbit that is just outside the stellar surface. 
The closest orbit to the star has been chosen because it is a 
characteristic orbit around the star. Another characteristic 
orbit that some stars have and belongs to the first category, 
is the orbit where the nodal precession becomes zero. This 
is a novel property that NS spacetimes have that rotating 
black hole (Kerr) spacetimes don’t exhibit. At that location 
the observables that we consider are the orbital frequency 
and the periastron precession frequency. 

For the second category of observables, we consider the 
accretion disc efficiency, the temperature distribution of an 
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accretion disc and in particular the maximum temperature, 
the photon energy where the spectrum has the maximum 
integrated luminosity, i.e., the photon energy where there is 
the maximum energy output from the disc, and finally the 
maximum integrated luminosity itself. The latter quantity 
turns out to be equivalent to the accretion efficiency, as one 
would expect. 

Finally we briefly discuss how could one use the vari¬ 
ous observables to perform consistency checks of the various 
proposed phenomenological models that relate astrophysical 
observables, such as QPOs, to properties of the geodesics 
around NSs, or how well accretion discs for example are 
described by a thin radiatively efficient disc model. Further¬ 
more we discuss how could one measure the first three mul¬ 
tipole moments of a NS by using observations of different 
observables. The latter application could be used to con¬ 
strain and maybe eventually determine the EoS of nuclear 
matter at supra-nuclear densities. 

We should emphasize again that the main underlying 
assumption in our analysis is that the motion of particles or 
fluid elements is geodesic. This implies some restrictions for 
the systems for which this analysis is applicable. One such 
restriction is that the central object should have a low mag¬ 
netic field, because in a different case, the large scale mag¬ 
netic field could modify the orbital motion of fluid elements 
in the inner part of the disc, that would then deviate from 
being geodesic. Another related restriction would be that the 
disc is in a relatively quiet state with a moderate accretion 
rate, so that the steady state thin accretion disc model is 
applicable. These requirements impose selection criteria for 
the LMXB systems that one could use. 

But one could go beyond the analysis presented here, by 
assuming a different model for accretion or different mecha¬ 
nisms for the QPOs. Even in that case, the main idea behind 
this analysis would remain the same, i.e., that one could use 
an analytic background spacetime that can be parameterized 
by only two parameters in order to express astrophysical ob¬ 
servables in terms of only these two spacetime parameters. 
The difference in that case though would be that in contrast 
to the relativistic precession model for QPOs and the radia¬ 
tively efficient thin disc model for the accretion, one would 
have to include some more physical parameters that would 
depend on the particular model. 

1.2 Plan of the paper 

The plan of the paper is as follows, in section 2 we briefly 
present the analytic spacetime that describes the exterior of 
NSs and discuss how we implement the universal relations 
between the moments to arrive to a universal NS space- 
time. In section 3 we first discuss the observables related to 
geodesic motion and then we discuss the observables related 
to the radiation emitted by an accretion disc around a NS. 
In section 4 we discuss how could one apply the results of 
section 3 to testing the various models for the astrophysical 
phenomena observed around NSs and how one could mea¬ 
sure the first three multipole moments of the NS. We end 
with the conclusions in section 5, while there is some further 
discussion in the appendix on an approximately universal re¬ 
lation for the equatorial radius of NSs (appendix A), on an 
approximately universal relation for the rotation frequency 
of NSs (appendix B), and on the properties of the nodal pre¬ 


cession frequency for NSs (appendix C). For the calculation 
we use geometric units, where G = c = 1 and the mass is 
given in km, unless some other unit is specified (usually for 
frequencies, temperature and energy). 


2 THE SPACETIME AROUND NSS 

Determining the spacetime around rotating compact fluid 
configurations is a difficult problem that can be solved an¬ 
alytically only in the slow rotation limit (although one still 
has to integrate numerically the interior of the star). If 
one were to tackle the problem of rapidly rotating con¬ 
figurations one would have to resort to a fully numeri¬ 
cal solution. The slow rotation solution is given by the 
Hartle & Thorne (1968) approach and there has been a 
lot of work extending this formalism up to fourth order 
in rotation (see work by Berti et al. (2005) and Yagi et al. 
(2014)). Similarly there has been a lot of work on nu¬ 
merical algorithms for solving the full Einstein field equa¬ 
tions for axisymmetric spacetimes around rotating fluid 
configurations (for example see Stergioulas & Friedman 
(1995) and for more details the review by Stergioulas 
(2003)). Additionally to these two approaches there have 
been investigations of analytic axisymmetric spacetimes 
that can match and accurately describe the exterior 
of rotating fluid configurations (see for example work 
by Stute & Camenzind (2002); Berti & Stergioulas (2004); 
Pappas (2009); Teichmiiller et al. (2011); Pachon et al. 
(2006); Pappas & Apostolatos (2013)). In particular, in 
Pappas (2009); Pappas & Apostolatos (2013) it was shown 
that the two-soliton analytic solution of Manko et al. (1995) 
can describe accurately the exterior spacetime of a NS for 
arbitrary rotation, rapid or slow. For our purposes here, we 
will use for the spacetime exterior to NSs the two-soliton 
analytic spacetime. 

The two-soliton solution is a four parameter solution 
of the vacuum Einstein’s field equations that is generated 
using the Ernst potential and the algorithm developed by 
Sibgatullin (1991); Manko & Sibgatullin (1993); Ruiz et al. 
(1995); Manko et al. (1995). 

The vacuum region of a stationary and axially symmet¬ 
ric spacetime in GR can be described by the line element 
first used by Papapetrou (1953) 

ds 2 = —f (dt — Lodcj)) 2 + [e 27 (dp 2 + dz 2 ) + p 2 d<(> 2 ] , ( 1 ) 

where /, uj, and 7 are functions of the Weyl-Papapetrou 
coordinates ( p,z ). By introducing the complex potential 
£{p,z ) = f(p,z) + ii/i(p,z), Ernst (1968) reformulated the 
Einstein field equations in the form of the complex equation 

(. Re(£))V 2 £ = V£ • VS, (2) 

where V and V 2 are respectively the gradient and the Lapla- 
cian in flat cylindrical coordinates (p,z,<j>). 

One can generate a solution of the Ernst equation by 
making a choice for the Ernst potential along the axis of 
symmetry of the spacetime, in the form of a rational function 

£{p = 0,z) = e(z) = (3) 

where P(z),R(z) are polynomials of z of order n with com¬ 
plex coefficients in general. 
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The vacuum two-soliton solution (proposed by 
Manko et al. (1995)) is a special case of the previous 
general axisymmetric solution that is obtained from the 
ansatz 


(z — M — ia)(z + ib) — k 
[z T M — ia)(z + ib) — fc ’ 


( 4 ) 


where all the parameters M, a, k, b are real (for more de¬ 
tails on the algorithm for generating solutions and for the 
two-soliton spacetime in particular see Manko et al. (1995); 
Pappas & Apostolatos (2013)). 

The first few multipole moments of the two-soliton 
spacetime are expressed in terms of the four parameters of 
the Ernst potential as, 


M 0 


M 4 


Ji 


M, M 2 = -(a 2 - k)M, 


a - (3a 2 - 2 ab + b 2 )k + k 2 - ifcM 2 
aM, Jo = — [a 3 — (2a — b)k]M, 


M 


( 5 ) 


where Mo = M is the mass, M 2 = Q is the quadrupole 
moment, Ma is the mass hexadecapole, Ji = J is the angular 
momentum, and J 3 is the spin octupole moment. 

It was shown by Pappas & Apostolatos (2013) that 
if one were to chose the parameters M, a, k,b in such a 
way so that the first four non-zero moments of the two- 
soliton spacetime were equal to the moments of the numer¬ 
ical spacetime (calculated using the prescriptions presented 
in Pappas & Apostolatos ( 2012 ) and further elaborated in 
Yagi et al. (2014)), then the analytic solution would be an 
accurate match for the numerical solution and that it would 
capture all the properties of the geodesics that are relevant 
to the astrophysical processes that take place in LMXBs. 

Furthermore, Pappas & Apostolatos (2014) have found 
that for NSs (and this has been extended to quark stars as 
well by Yagi et al. (2014)) the first multipole moments are 
not all independent between them. In particular one can ex¬ 
press the first moments that are higher than the quadrupole 
in terms of the quadrupole. Specifically, if we define the re¬ 
duced moments as 


M-n 


M n 

jn-M n + 1 ’ 


Jn 


J n 

j n M n+1 ’ 


( 6 ) 


where j is the spin parameter defined as J/M 2 , then the 
spin octupole and the mass hexadecapole of a NS will be 
related to the quadrupole by relations of the form 


y = A + Bax " 1 +B 2 x V2 , (7) 

where y can be either y/— J 3 or yj Ma and x is V-m 2 . 
Therefore the first higher moments of a NS spacetime can 
be expressed in terms of only three parameters, the mass M, 
the angular momentum J, and the quadrupole M 2 = Q. 

Additionally, as it was shown by Pappas & Apostolatos 
(2014) and Yagi et al. (2014), these relations between the 
moments are independent of the EoS which means that if one 
were to use these expressions to produce a two-soliton space- 
time, then one could have an EoS independent description of 
the spacetime around a NS. To clarify the last statement, us¬ 
ing the universal relations between the moments, one could 
construct a spacetime metric parameterized in such a way, 
so as to be suitable to describe any NS having any one of 
the realistic EoSs that we have at our disposal. This EoS- 
independent (or -universal) description of the NS spacetime 
is what we aim to have here. 


Since for the construction of the spacetime we have only 
four parameters which are fixed by the first four moments, 
M, J, M 2 , and J 3 , we will only need the relation between the 
quadrupole and J 3 given by Pappas & Apostolatos (2014), 
i.e., 

y = —0.36 + 1.48 a: 0 ' 65 , ( 8 ) 

where y = y/— J 3 and x = y/—M 2 . For the description of the 
spacetime and the various properties of the geodesics that 
we will calculate, we will use as a first parameter the dimen¬ 
sional mass M, expressed in units of km, and two additional 
dimensionless parameters which will be the spin parameter 
j = J/M 2 and the reduced quadrupole a = —M 2 /(j 2 M 3 ). 


3 ASTROPHYSICAL PROPERTIES AND 
OBSERVABLES 

As it was elaborated in the previous section, we now have a 
three parameter description of the spacetime around a NS, 
which is independent of a particular choice for the EoS of 
the matter inside the NS. The three parameters that enter 
the description are the mass M that gives the scale to all the 
quantities and the two dimensionless parameters j, which is 
the spin parameter and is a measure of the rotation of the 
NS, and a, which is the reduced quadrupole and is a mea¬ 
sure of the deformability due to rotation. For this reason we 
will also call the parameter a the quadrupolar deformability 
parameter. Since the mass M is a scale of the spacetime and 
its properties, by appropriately reducing or normalizing the 
various quantities we can end up with a two-parameter de¬ 
scription of all the astrophysical observables, as we will see 
in what follows. 


3.1 Equatorial geodesics in an axisymmetric 
spacetime 

We start with a brief discussion of geodesics in an axisym¬ 
metric spacetime. An axisymmetric spacetime given by the 
line element in eq (1) admits two Killing fields which de¬ 
fine the symmetries with respect to translations in time and 
rotations with respect to an axis. One can take advantage 
of these symmetries and the corresponding integrals of mo¬ 
tion to study geodesics in the spacetime. The first integral 
of motion is the energy E, given as 

E =-paC = ~Pt = y(^-gu^ - , (9) 

where t is the coordinate time and r is proper time. The 
second integral of motion is the angular momentum L in 
the direction of the symmetry axis, given as 

L = PaV a = p<)> = p (gtt/j// + (40) 

In addition to the two previous equations we have from the 
normalization of the four-velocity the equation 


- 1 = 


9tt \/t) +29t *(£) (d/ ]+i 


+9pp l d-7 


, dz 
+ 9zz dJ 


( 11 ) 
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In this equation, one can define the angular velocity fi = 
Then for the circular and equatorial orbits, eq. (11) defines 
the redshift factor between coordinate and proper time, 

^ ~ gtt , (12) 



and the energy and the angular momentum take the form, 


e- e - 

gtt gtci>Q 

(13) 

g 

yj gtt 2igt(f)Q gcfxfrQ^ 

l- l - 

gt* + g 

(14) 

g 

y/ gtt 2 gt(f)Q 


where we have introduced the energy and angular momen- 

, , 2 

turn per unit mass. From the conditions, = 0, = 

0, z = 0 and 4| = 0 for circular equatorial orbits, and 
the equations of motion obtained assuming the Lagrangian, 
C = \g a bx a x b , the orbital angular velocity can be calculated 
to be, 


9t<t>,p + y/(9t<t>,p) 2 9tt,/>gp 

9<t>4'iP 


This is the orbital frequency of a particle in a circular orbit 
on the equatorial plane and the commas indicate partial 
derivatives with respect to the corresponding coordinates. 

Equation (11) can take a more general form in terms of 
the constants of motion, 

( dp\ ( dz\ _ E 2 g^rj, + 2ELgtcj > + Li 2 gtt 

9pp \dr) 9zz \dr) {9 h) 2 ~ gttg^ 

= Veff. (16) 

With equation (16) we can study perturbations around cir¬ 
cular equatorial orbits. If we assume small deviations of the 
form, p = p c + 5p and 2 = Sz, then we obtain the perturbed 
form of equation (16), 


9pp 




1 d 2 Veff 

2 dp 2 


(<5p ) 2 


, 1 d 2 Veff 

2 dz 2 


(Sz) 2 (17) 


This equation describes two harmonic oscillators with fre¬ 
quencies, 


- 2= g pp d 2 v eff 

p 2 dp 2 
2 _ g” 8 2 V eff 
z ~ 2 dz 2 


(18) 

(19) 


where everything is evaluated at the corresponding equato¬ 
rial circular orbit. The first frequency is the frequency of ra¬ 
dial oscillations around the radius of the circular orbit, while 
the second frequency is the oscillation over and under the 
equatorial plane. The differences of these frequencies (cor¬ 
rected with the redshift factor (12), i.e., K a = (dr / dt)n a ) 
from the orbital frequency, fl a = Q, — K a , define the pre¬ 
cession frequencies, where the oscillation frequencies K, a are 



Figure 1. The figure shows contours of constant Risco/M on 
the parameter plane of ( j,a ). The range of the parameters (j, a) 
has been chosen so as to include all physically relevant NS models, 
i.e., models with masses in the range between the maximum stable 
mass and a minimum of around lMg with rotation from zero up 
to the maximum allowed limit. 


given in terms of the metric functions as, 

\ " / ,aa 

— 2 {g t t + gt(/>£l){gt(i> + ( 20 ) 

V P / ,aa 

+ ( flt ^ + S 4><#>^) 2 

where the index a takes either the value p or 2 to express the 
frequency of the radial or the vertical perturbation respec¬ 
tively and the expressions are evaluated on the equatorial 
plane 2 = 0. The position where n 2 becomes zero is the 
location of the ISCO. 

In what follows we are going to look into corotating 
orbits and all the relevant discussion will be on this type of 
orbits. 



3.2 Location of the ISCO 

The first property of the geodesics around a NS that we will 
examine is the location of the ISCO. The ISCO is located at 
the position where k 2 p — 0. Assuming the two-soliton ana¬ 
lytic spacetime with the moments chosen as it was described 
in section 2, utilizing the universal relations between the 
moments, we can have by solving k 2 = 0 the location of 
the ISCO for a general NS spacetime independently of any 
specific EoS. 

Using as free parameters the spin parameter j and the 
quadrupolar deformability parameter a, we can calculate the 
reduced circumferential radius of the ISCO, Risco/M. Fig. 
1 shows the contours of constant Risco /M. Something that 
we should point out here is that the contour of Risco/M = 
6 , as the spin parameter j decreases, goes up asymptotically 
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Figure 2. The figure shows contours, on the parameter plane of 
(j, a), of constant relative difference between the Risco of the 
spacetime and the Rjgco calculated with the Shibata &: Sasaki 
formula up to third order in the expansion, i.e., up to the order 
that includes the spin octupole. 

to the axis of a and when j = 0 it spans the entire axis of 

a. 

One thing that would be interesting to explore here, 
since we have the EoS-independent ISCO, is the generic ac¬ 
curacy of the formula developed by Shibata & Sasaki (1998) 
that gives the radius of the ISCO as an expansion on the mul- 
tipole moments. The formula developed by Shibata & Sasaki 
gives the ISCO radius as an expansion in the multipole mo¬ 
ments of the spacetime going up to the mass hexadecapole, 
M 4 . Since we are using here expressions that are taking con¬ 
sistently into account the moments up to the spin octupole, 
J 3 , and since Berti & Stergioulas (2004) have demonstrated 
that it is better to be consistent in the order of the expan¬ 
sion that is used in the Shibata & Sasaki formula, we will 
use the formula up to the third order in the expansion, i.e., 

RTsco = 6 Af(l - 0.54433j - 0.22619/ + 0.17989a/ 
-0.23002/ + 0.26296a/ - 0.05317/0),(21) 

where a is the reduced quadrupole in our formalism and 
/3 = —Ja/(/M 4 ) is the reduced spin octupole. The pa¬ 
rameter /3 in this formula will be given in terms of the 
parameter a, as it was discussed at the end of section 
2, having therefore the ISCO in our two dimensional pa¬ 
rameter space, (j, a). In Fig. 2 we show a contour plot of 
the relative difference between the Risco and the Risco > 
A Risco = ( Risco — Risco)/Risco- As one can see from 
the plot, the Shibata & Sasaki formula is very accurate for 
small values of the spin parameter j and it remains accurate 
for a large part of the parameter space. The relative differ¬ 
ence goes over 10 per cent only on the upper right corner 
of the parameter space, where we have the rapidly rotating 
NSs that have also a relatively large deformability. This is 
the behaviour that one would expect, since for these models 
the deformation is large enough to make the higher order 
moments more relevant for the calculation of the ISCO. 


Figure 3. The figure shows contours, on the parameter plane of 
( j,a ), of constant relative difference between the Risco of the 
spacetime and the equatorial circumferential radius of the star 
Req, i.e., ( Risco — Req)/Risco ■ The contours with a positive 
value correspond to a region of the parameter space where the 
stars have their surfaces under the ISCO. The negative values 
correspond to stars that have engulfed their ISCO. 

One last thing that we should consider at this point is 
for which part of the parameter space, the position of the 
ISCO is outside of the surface of the star, i.e., for which 
region of the parameter space the equatorial circumferential 
radius of the stellar surface R eq is smaller than the Risco- 
In the appendix A we discuss how one could have an EoS- 
independent description of the equatorial radius of a NS in 
terms of the two dimensional parameter space that we are 
using. Fig. A1 gives the reduced circumferential equatorial 
radius of a NS in terms of the parameter space (j , a). What 
we would like to see is for which part of the parameter space 
the inequality R eq < Risco holds. 

For that reason we plot in Fig. 3 the relative differ¬ 
ence between the Risco of the spacetime and the equa¬ 
torial circumferential radius of the star R eq , i.e., ( Risco — 
Req)/Risco- Positive values correspond to models that have 
their ISCO outside the stellar surface, while negative values 
correspond to models that have their ISCO buried under the 
surface of the star. Fig. 3 shows that when the deformabil¬ 
ity of a star is large enough, then the ISCO is buried under 
the surface. But as we can see by comparing Figs 3 or A1 
against Fig. 1, the radius of the star is not that much larger 
than the ISCO and it can be up to around 50 per cent of 
the ISCO, i.e., for a star that has its ISCO in the range of 
6 — 7 M the surface will be in the range of 10 — VIM. 

This is something that one needs to take into account 
when considering observables from orbiting matter around 
a compact object. The presence of the surface will constrain 
the observed orbital and precession frequencies that a fluid 
element can have as it orbits around a NS and it will also 
constrain the appearance of an accretion disc that will either 
terminate at the ISCO or at the surface, depending on the 
model. 

One final point that one could make here is that, from 
comparing Figs 2 and 3, one can conclude that for the part 
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Figure 4. The figure shows contours, on the parameter plane 
of (j, a), of constant reduced frequency M X vi SCO = M X 
^ISCO /(2tt) in units of km X kHz. 

of the parameter space that an ISCO exists outside the NS, 
the Shibata & Sasaki formula is extremely accurate (better 
than 1 per cent). 

3.3 Orbital frequency flisco at the ISCO 

The second property of the geodesics that can be a po¬ 
tential astrophysical observable is the orbital frequency at 
the ISCO, SI/sco, or rather the corresponding frequency 
v = Q,/(2n). One would expect that this frequency might 
appear in QPOs associated to the orbital motion of fluid 
elements in the accretion disc. In particular it might have 
distinctive characteristics, since at that point the accreting 
material plunges to the compact object and the disc termi¬ 
nates. In Fig. 4 we have plotted contours of the normalized 
frequency M x visco in units of km x kHz. 

In the previous section though, we talked about the lo¬ 
cation of the ISCO with respect to the surface, and we saw 
that for some models it is possible that the surface of the 
compact object overcomes the location of the ISCO. This 
results to having the available space for orbital motion ter¬ 
minate before the ISCO, at the radius of the surface. This 
means that since we have to take also into account the sur¬ 
face, the figure for the frequency will be modified. We can 
therefore plot essentially the maximum orbital frequency, 
i/^ ax , which is the orbital frequency of the closest stable cir¬ 
cular orbit to the surface of the compact object. The relevant 
contour plot is given in Fig. 5, where we also indicate the 
models for which the ISCO radius is equal to the surface ra¬ 
dius, Risco = Req ■ Below that curve the orbital frequency 
plotted is that of the ISCO, while over that curve the or¬ 
bital frequency plotted corresponds to the orbital frequency 
at the surface of the star. 

It is tempting at this point to create a combined figure 
of the contours for the maximum orbital frequency and the 
rotation frequency of a NS, presented in the appendix B. 
Both of these frequencies are potentially observable and are 
probably the most prominent observables one might have 



Figure 5. The figure shows contours, on the parameter plane of 
(j , a ), of constant reduced frequency M X i/^ ax in units of km X 
kHz The dashed line that intersects the solid lines represents the 
models for which the ISCO radius is equal to the surface radius. 
Below that curve the orbital frequency plotted is that of the ISCO 
which is outside the surface of the star. Above that curve the 
orbital frequency plotted corresponds to the orbital frequency at 
the surface of the star, since the ISCO is under the surface. 


from an accreting rotating NS. We present this figure at 
this point because it displays an interesting physical charac¬ 
teristic of rotating NSs which this parameterization nicely 
demonstrates. 

Fig. 6 shows the contours of constant reduced maximum 
frequency M x overlaid against the contours of con¬ 

stant rotation frequency of NSs. The interesting feature of 
this figure is the vertical dashed curve that is located in the 
region above the curve of Risco = Req, he., in the region 
where the dotted curves of maximum orbital frequency cor¬ 
respond to the orbital frequency at the radius of the surface. 
One notices that along the vertical dashed curve, curves of 
constant NS rotation frequency intersect curves of the same 
constant orbital frequency. This signifies the models that 
are rotating at the Kepler limit, i.e., at a frequency such 
that a fluid element on the surface of that star would have 
to rotate at a Keplerian orbital frequency. This means that 
the fluid elements at the surface of such a configuration are 
not gravitationally bound to the star and all models to the 
right of the vertical dashed curve, i.e., models that rotate 
more rapidly than the Keplerian orbital frequency, can not 
be physically realized. The vertical dashed curve therefore 
represents a physical boundary to our parameter space for 
NSs, the boundary of maximally rotating models. It is in¬ 
teresting to note that the Kepler limit is at j ~ 0.72 and is 
quite independent of the parameter a. We should also note 
that since our description applies for all EoSs, the limit that 
we get for the spin parameter is EoS-universal. This result is 
reminiscent of the universal empirical formula for the Keple¬ 
rian frequency by Lattimer & Prakash (2004); Haensel et al. 
(2009), but it is not clear how one could get from that for¬ 
mula to this result. 
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Figure 6. The figure shows the contours, on the parameter plane 
of (j, a), of constant reduced frequency M X i/j£ ax (dotted curves) 
overlaid against the contours of the rotation frequency of NSs 
(solid curves), in units of km X kHz. The figure shows also the 
dashed curve of Risco = Req and another vertical dashed curve 
that represents the models that are rotating at the Kepler limit. 
This becomes evident if one notices, that the curves of constant 
stellar rotational frequency intersect the curves of constant orbital 
frequency at the surface, that have the same value. The region to 
the right of that curve corresponds to physically not realisable 
stellar models. 

3.4 Vertical oscillation frequency kisco at the 
ISCO 

Another geodesic property that we could associate with the 
location of the ISCO is the nodal precession of a slightly off- 
equatorial orbit, v z = (fi — K z )/(2n) at the location of the 
ISCO. The behaviour of the nodal precession frequency of 
NSs is quite different from the behaviour of the correspond¬ 
ing Kerr BHs. In particular the nodal precession exhibits 
some very interesting features that can not be found in ro¬ 
tating BHs and essentially demonstrates the fact that the 
exterior spacetime of NSs deviates from that of Kerr both 
quantitatively and qualitatively. To demonstrate this, we are 
plotting in Fig. 7 the normalized nodal precession frequency, 
M x v z at the location of the ISCO. 

The feature that stands out is the fact that the nodal 
precession is positive at the ISCO for some models while 
for other models it is negative. This effect, which had been 
noticed by Pappas (2012) in connection to QPO observa¬ 
tions and NS properties, is further discussed in appendix 
C. This means that if there were a slightly inclined accre¬ 
tion disc around a NS that demonstrated negative preces¬ 
sion frequency at the ISCO, then that accretion disc would 
have the inner region precessing in one direction while the 
outer region would precess in the opposite direction (see 
discussion in appendix C). This effect could be of signifi¬ 
cance to the modelling of accretion discs and could have even 
more spectacular consequences than the effects explored by 
Nixon et al. (2012) for accreting BHs. 

In this case as well as in the previous cases, one needs 
to take into account the effect that the surface of the NS 
will have, since the presence of the surface in some cases 


Figure 7. The figure shows contours, on the parameter plane 
of ( j,a ), of constant reduced nodal precession frequency M X 
(Q z /2 tc)isco i n units of km X Hz. 



Figure 8. Same plot as in Fig. 5, but for the nodal precession 
frequency. Again the horizontal dashed line represents the models 
for which Risco = Req, while the vertical dashed line represents 
the models at the Kepler limit. 


will not allow matter to access radii all the way down to the 
ISCO. Therefore, if we focus only on accessible precession 
frequencies the previous picture will be modified. In Fig. 8 
we show the contours of constant nodal precession frequency 
with the modification that above the horizontal dashed line 
of Risco = Req, the frequency corresponds to the location 
of the surface. We can see from Fig. 8 that even though the 
presence of the surface changes things, there are still models 
that can have a negative precession frequency at the inner 
part of the spacetime close to the surface of the star. 
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Figure 9. Same plot as in Fig. 5, but for the orbital frequency 
(black long dashed curve) and the periastron precession frequency 
(red dash-dotted curve) at the location where the nodal preces¬ 
sion becomes zero. The frequencies are in units of km x kHz. The 
almost diagonal dotted line indicates the models for which the 
radius of the surface coincides with the radius where the nodal 
precession is zero, while the vertical dashed line represents the 
models at the Kepler limit. Also, the vertical axis has been mod¬ 
ified in order to include only the relevant region of the figure. 
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3.4-1 Region of zero nodal precession 

The effect discussed in the previous section, i.e., that the 
nodal precession near the surface of the star is negative, 
means that there is a radius further out from the surface 
where the nodal precession is zero. In the case that we have 
an accretion disc, that region is a place where there might be 
something significant taking place, such as the breaking of 
the disc for example. Therefore it would be of interest to see 
what is the behaviour at that point of the other two frequen¬ 
cies, i.e., the orbital frequency and the periastron precession 
frequency. Contours of these two frequencies are plotted in 
Fig. 9. The dotted diagonal line indicates the models for 
which the radius of the surface coincides with the radius 
where the nodal precession is zero. Below that curve there 
is no zero nodal precession outside the star. 


3.5 Thin accretion disc 

The relativistic model of a radiatively efficient thin accre¬ 
tion disc was developed by Novikov & Thorne (1973) and 
Page & Thorne (1974). The disc is thin in the sense that 
the vertical size of the disc, H is much smaller than the 
characteristic radial length of the disc. It is also radiatively 
efficient in the sense that all the gravitational energy lost 
by the accreting matter due to viscous stresses is efficiently 
radiated away or transported outwards, therefore there is 
no heating of the disc. The disc is also considered to be 
in a hydrodynamic equilibrium and in a steady state where 
the accretion rate Mo is constant. Furthermore, the pressure 
gradient and the heat flow within the disc are negligible in 
the radial direction, while heat can flow in the vertical di¬ 


rection (for details on the assumptions see Page & Thorne 
(1974)). 

This construction results in an accretion disc where the 
plasma moves in equatorial Keplerian orbits with angular 
velocity 12, energy per unit mass E and angular momentum 
per unit mass L, slowly inspiraling towards the central ob¬ 
ject releasing gravitational energy that is radiated from the 
two faces of the disc. In the end, the structure equations 
of the disc are given by the mass conservation, the energy 
conservation and the angular momentum conservation. In 
this construction, everything is integrated along the thick¬ 
ness of the disc and the metric is assumed to be on the 
equatorial plane z = 0, while deviations from the equatorial 
plane behave as (z/p) 2 , due to equatorial symmetry, and are 
considered negligible. 

The energy momentum tensor for the accreting matter 
fluid has the form, 

T= p 0 u p u + 2u ( V ) + t^ (22) 

where u p is the four-velocity of particles, q p is the energy 
flow four-vector, po is the rest mass density, and is 
the stress tensor (the specific internal energy II is negligi¬ 
ble according to the initial assumptions). The first struc¬ 
ture equation comes from the conservation of the rest mass, 
V M (po« M ) = 0, and it essentially states the fact that the 
time-averaged mass accretion rate is independent of the ra¬ 
dius, i.e., 

Mo = —2-iis/—gEu p = const., (23) 

where g is the determinant of the metric (in our case the 
metric (1), which is calculated on the equatorial plane), u p 
is the radial component of the four-velocity, and E is the 
surface density defined as, E = f H H {po)dz, where (po) is the 
rest mass density averaged in time and over the azimuthal 
angle 4>. 

The second equation of structure comes from the angu¬ 
lar momentum conservation and is, 

[MoL - 2ny/=jW$] = 4tt ^FL, (24) 

where F is the radiated energy flux from one of the surfaces 
of the disc and W p is the time-averaged torque per unit 
circumference at a particular radius, due to the stresses. The 
latter is defined as W£(p) = J H H {t^,)dz. The last equation 
of structure comes from the energy conservation and is, 

[m 0 .E - 2-K^^nW^ = 4t t^/^FE. (25) 

To these equations we need to add the fundamental relation 
for circular geodesics, i.e., 

E, p = QL, P . (26) 

By making the redefinitions 

/ = An^F/M 0 , w = 2n^gW p /M 0 , (27) 

the angular momentum and energy conservation equations 
can take the form 


C L-w), p = fL , (. E-Shv), p = fE. (28) 

From these two equations one can obtain the algebraic rela¬ 
tion for w, 


w - 


E-nL 


(29) 
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and by substituting this, in the previous equations, we reach 
after some algebra to the integrable equation for /, 


( (E~9L) 2 \ 

V (-n,„) J ), p 


= (E-9L)L, P . 


(30) 


For this equation to be integrated, one needs to choose the 
appropriate boundary conditions. The usual assumption is 
that the disc terminates at the ISCO, where matter plunges 
towards the central object. This means that at the ISCO, 
which is the inner edge of the disc, one assumes a zero torque 
condition. With this boundary condition the resulting equa¬ 
tion for the flux is 


/ 


n, P 

{E - 9L) 2 


f 

J PISCO 


(E-nL)L y dp. 


(31) 


As we have seen from our discussion so far, for NSs this 
is not the only possibility. For some models the surface of 
the star overcomes the location of the ISCO, therefore in 
these cases the disc should terminate on the NS surface. 
Apart from choosing the inner edge of the disc to be at the 
radius of the surface, we also need to choose and appropri¬ 
ate boundary condition for the torque. As it is discussed 
by Zimmerman et al. (2005), in the region where the accre¬ 
tion disc touches the surface, a boundary layer forms that 
causes the orbital frequency of the accreting material to de¬ 
crease until it reaches that of the surface of the rotating star. 
In that region, we have a maximum of the orbital frequency 
where there are no viscous shear stresses and therefore there 
is zero torque. The orbital frequency at the maximum is es¬ 
sentially the maximum geodesic orbital frequency, i.e., the 
orbital frequency at the radius of the surface. In this case 
therefore, the result for the flux is 


/ = - 755 


n„ 


(. E-9L ) 2 


/' 

J P.s 


(E - 9L)L y dp , 


(32) 


where the subscript s indicates the values at the surface of 
the star. 

Finally, the emitted flux from the surface of the accre¬ 
tion disc as a function of the radius will be given by the 
expression 


np) = (33) 

Since the disc is considered to be in thermal equilibrium, the 
surface of the disc at each radius radiates as a blackbody. As 
a consequence, the emitted flux is related to the blackbody 
temperature by the Stefan-Boltzmann law 

F(p) = <jT 4 (p), (34) 

where a is the Stefan-Boltzmann constant. 3 Therefore from 
the distribution of the radiated flux we obtain the tempera¬ 
ture distribution of the disc. 

We should note here how the different quantities scale 
with the mass of the central object, since we will need it for 
our analysis. As we have already seen, the radial distances 
scale proportionally to the mass, p = Mp, while the frequen¬ 
cies scale as 9 = M~ 1 9. Similarly, the orbital energy and 
angular momentum per unit mass scale as E = M°E and 


3 In the relevant units for us, the constant is a = 1.562 X 
10 —54 km“ 2 K~ 4 



Figure 10. Same plot as in Fig. 5, but for the efficiency p in 
percentile. The efficiency is independent of the mass of the cen¬ 
tral object and it is evaluated either at the radius of the ISCO 
or at the radius of the surface of the star. The almost horizon¬ 
tal dashed line indicates the models for which the radius of the 
surface coincides with the radius of the ISCO, while the vertical 
dashed line represents the models at the Kepler limit. 


L = ML. The square root of the determinant of the met¬ 
ric, which is essentially a radius (the metric (1) is given in 
cylindrical coordinates), also scales linearly with the mass. 
Therefore, the redefined flux / will scale as f(p) = M~ 4 f(p), 
while the flux o ver the mas s accretion rate will scale as 
F(p)/Mo = M~ 2 ^F(p)/Moj ■ From these scalings we have 


that, modulo the accretion rate, the temperature scales as 
T(p) = M~ 1/2 T(p). Actually it would be convenient to add 
the accretion rate to the scaling, in the sense that we could 
express the Stefan-Boltzmann law as F/Mo = crT 4 , where 
T = (Mo)~ 1,/4 T is an accretion rate reduced temperature, 
and get the scaling, T = M~ 1 ^ 2 (Mo) 1 ^ 4 T. 


3.6 Emitted efficiency p of accretion 

A first measure of the luminosity of the disc can be given by 
the efficiency p, defined as 


P = 1-e\ , (35) 

I in 

where E\ is the energy per unit mass of the circular orbit at 

I in 

the inner edge of the accretion disc, which can be either the 
location of the ISCO or the surface of the star, as we have 
seen so far. The efficiency p is essentially a measure of the 
energy that has been radiated away to infinity from the disc 
surface per unit mass of infalling material, assuming that 
all the photons reach infinity. Since E scales independently 
of the mass of the central object, the efficiency will also be 
independent of the mass. 

In Fig. 10 we present the contours of constant efficiency 
p, evaluated at the inner edge of the disc. The figure is sep¬ 
arated in two parts, the lower part under the dashed curve 
of Risco = Req, where the inner edge of the disc is at the 
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ISCO and the upper part where the inner edge of the disc 
is at the surface of the star. The values are given per cent. 
One can get the luminosity that corresponds to the given 
efficiencies by multiplying the efficiency to an accretion rate 
Mo- 


3.7 Maximum disc temperature 


In the context of the thin disc model, one can calculate the 
emitted flux profile and from that the temperature distribu¬ 
tion for a given model, as we have previously discussed. A 
quantity that would be characteristic of the emitted spec¬ 
trum can be considered to be the effective temperature that 
corresponds to the maximum flux. Therefore we will present 
here the behaviour of the maximum effective temperature of 
thin accretion discs around NSs. 

We could present the results for the fluxes and the ef¬ 
fective temperature distributions in terms of the mass and 
accretion rate reduced quantities that we discussed in the 
end of section 3.5, but instead we will prefer to do the cal¬ 
culations with respect to a fiducial accretion rate, which we 
will chose to be safely below the Eddington limit and close 
to a plausible outflow value from a companion star in a 
LMXB. Specifically we will choose Mo = 10~ 12 Mgyear _1 
(see Shakura & Sunyaev (1973)). The results in this way 
will be more intuitive. Apart from that we will maintain the 
mass-reduced quantities. Furthermore, for reasons of more 
compact and clearer ilustration, in some cases instead of 
giving the temperature in terms of degrees Kelvin, we will 
prefer to give it in terms of electron volts by assuming the 
correspondence of leV/fes = 11, 604K (fes is the Boltzmann 
constant). The resulting temperature will be in eV km 1 / 2 , 
since the temperatures that we are calculating are mass- 
reduced. In any case, one can express the temperatures pre¬ 
sented here, in terms of any accretion rate and mass of the 
central object using the following relation 


T 

eV 


= 1.477 


M_ 

M© 

T 


- 1/2 


Mo 


lO _ 12 M 0 year _1 


leV km 1//2 


1/4 


(36) 


This expression can be interpreted to mean that the tem¬ 
peratures presented here are calculated having as a unit of 
accretion rate the rate of lO~ 12 M 0 year _1 and as a unit for 
the mass of the central object, either the geometric unit of 
lkm or the unit of ( 1 . 477 ) - 1 M 0 . 

As it was shown in section 3.5, for a given choice of the 
spin parameter j and the quadrupolar deformability parame¬ 
ter a, which specify a particular geometry for the spacetime, 
one can calculate the emitted flux of a thin and radiatively 
efficient disc around a NS as a function of the radius. Then 
from the Stefan-Boltzmann law (34) one can obtain the tem¬ 
perature profile of the disc as a function of the radius. Fig. 
11 shows such temperature profiles for different choices of a 
and j. In particular, assuming the aforementioned accretion 
rate, the plots show the mass-reduced effective temperature 
T = M 1 / 2 T as a function of the reduced circumferential 
radius R/M. The two top panels in Fig. 11 correspond to 
models that have an ISCO, while the two bottom panels cor¬ 
respond to models without an ISCO. What is apparent from 
comparing the top and the bottom panels is that the shape 


of the temperature distribution is affected by the location of 
the inner edge of the disc. In the two upper panels the inner 
edge is determined by the location of the ISCO and therefore 
it depends in a specific way on the rotation (changes mono- 
tonically with the spin), while in the two bottom panels the 
inner edge is determined by the location of the surface of the 
star which changes the behaviour with respect to rotation. 
The observed behaviour in the latter case becomes clear if 
one observes Fig. Al, where one can see that for a given 
value of a and for increasing spin parameter j, at first the 
radius of the star decreases and after a point it starts to 
increase again. 

A characteristic property of these distributions that 
could be of interest is the maximum effective temperature 
T^ ax . One could construct contour plots of that temperature 
in the parameter space of (j, a ). The contours for the maxi¬ 
mum effective temperature are given in Fig. 12. The plot in 
Fig. 12 also shows the curve that separates the models that 
have the ISCO outside the surface of the star from those 
that the surface has overcome the ISCO, as well as the limit 
of maximum rotation of the NSs, as previously discussed. 

The surface that gives the maximum effective tempera¬ 
ture as a function of the spin parameter j and the quadrupo¬ 
lar deformability parameter a, can be quite accurately de¬ 
scribed by a function of the form, 


T e ff ax = Bo + Boj + B 2 j 2 + (A 0 + Aij + A 2 j 2 )a ni 

+(Co + C 1 j+C 2 j 2 )a n2 , (37) 

where A 0 = -3.95445, Ai = -12.97699, A 2 = 35.90089, 
B 0 = 76.77519, B x = -44.52375, B 2 = 358.23145, C 0 = 
18.142198, Ci = 89.3791, C 2 = -322.8577, m = 1.107, and 
n 2 = 0.4425. In particular, the relative difference between 
the fitting function and the actual values is always lower 
than 2 per cent as one can see in Fig. 13, where we have 
plotted the relative difference in percentile. 


3.8 Emitted spectrum 

For a given temperature distribution of the accretion disc 
and the corresponding geometry of the spacetime around the 
central object, one could calculate the spectrum that a dis¬ 
tant observer would measure. A detailed calculation would 
have to involve various effects, such as that of the propa¬ 
gation of photons in the background spacetime, i.e., light 
bending, and other effects. One can find a detailed analy¬ 
sis for example in the work by Li et al. (2005). Here, for 
illustrative purposes, we will follow a less exhaustive calcu¬ 
lation, where we will assume that the disc is almost face on 
(ignoring inclination effects and light-bending) and we will 
also neglect non-thermal effects and limb-darkening (neglect 
colour correction and assume isotropic emission). Following 
this approximation, one can straightforwardly calculate the 
emitted spectrum of the accretion disc. 

As we have discussed in section 3.5, the accretion disc 
has a temperature distribution that is given as a function 
of the radius and the material is assumed to be in local 
thermodynamic equilibrium and thus the emitted radiation 
from each location on the surface of the disc is a blackbody 
radiation. That means that the local specific intensity of the 
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Figure 11. Indicative temperature distributions of the accretion discs for different values of the quadrupolar deformation a and the 
spin parameter j. The quantity plotted is the mass-reduced temperature T = M 1 ' 2 T as a function of the logarithm of the reduced 
circumferential radius R/M. One can see that for the values of a and j that correspond to models where the ISCO location is outside the 
stelar surface (two top panels), the temperature distributions are ordered with rotation. That is because the shape of the distribution and 
the maximum effective temperature are mainly determined by the location of the inner edge of the disc, which is given by the location 
of the ISCO, which in turn depends on the rotation. On the other hand, for the models that have values of a and j that correspond to 
models where the ISCO location is inside the star (two bottom panels), the spin of the NS is not the main factor that determines the 
shape of the distribution. The main factor in these cases is the location of the surface, which determines the location of the inner edge 
of the disc. These temperature profiles have been constructed under the assumption of an accretion rate of Mo = 10 — 12 M 0 year — 1 as 
discussed in the main text. 


emitted radiation will be 


Ie{Ve) = 


2hvi 


exp 


/ hv e \ 
l fesTeff J 



(38) 


where u e is the emitted photon frequency and h is the Planck 
constant. For an observer that is very far from the disc, the 
observed specific flux will be 

F„ ob3 = / -fofos(^obs)dfl 0 bs, (39) 

where I obs is the observed specific intensity, v ob s is the ob¬ 
served photon frequency and d0. o b s is the observed differen¬ 
tial solid angle. The quantities at the observer’s frame can be 
calculated in terms of the quantities on the emission frame 
on the accretion disc, taking advantage of the relativistic 
invariant (7„/i/ 3 ). The specific flux will be then, 


F^obs 



Vobs)d£l 0 b Sl 


(40) 


where Q = is a redshift factor. Instead of calculat¬ 

ing the specific flux (which depends on the distance of the 
source) it would be more convenient to calculate the corre¬ 
sponding isotropic luminosity 


L v = 4:TtD 2 F„, (41) 

where D is the distance of the source from the observer, 
which in the end is independent of the distance since dQ 0 t, s = 
dE/D 2 , where dE is the surface at the source. The final 
expression for the luminosity in terms of the observed photon 
energy (neglecting light-bending effects) will be, 


L 


£ obs 


8-7T cos i 


/ 


£ 3 obs y/=gdpd(j> 

[hc? h> Usfc) -'] 


(42) 


where i is the inclination angle, which is zero for a face-on 
disc and is assumed to be small in our approximation. Since 
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Figure 12. Same plot as in Fig. 5, but for the maximum mass- 
reduced effective temperature T^ ax in eV km 1 / 2 , under the as¬ 
sumption of an accretion rate of Mq = ! 0 12 Mg year - 1 as dis¬ 
cussed in the main text. 



Figure 13. Contours of relative difference per cent of the maxi¬ 
mum effective temperature with respect to the value given by the 
fitting function. The absolute value of the relative difference is 
always smaller than 2 per cent. 

we are neglecting light-bending effects, the redshift factor 
will be (see for example Bambi (2012)) 


g _ \/ 9tt 2g^fi 2 
1 + tip sin <j> sin i 

where Q is the orbital frequency of the emitting fluid element 
on the accretion disc. Therefore, to calculate the spectrum 
we need to integrate from the inner edge of the disc up to 
the outer edge. For our calculation, the inner edge will be 
either the location of the ISCO or the surface, depending on 
the model, as it was discussed in the previous sections. For 
the outer edge we will assume a radius of around 10 3 M. 


At this point we should make some comments on how 
the various quantities scale with the mass of the central ob¬ 
ject and the accretion rate. As we discussed in section 3.5, 
the effective temperature scales with the mass and the ac¬ 
cretion rate as T = M _1 / 2 (Mo) 1 / 4 T. This means that the 
energy of the emitted and observed photons will also scale 
as Shu = M~ 1 t 2 (M 0 ) 1 t 4 ‘e hv . If we introduce that to the 
expression for the luminosity we can see that it scales as 
Lu = M 1 / 2 (Mo) 3 / 4 L y , while the integrated luminosity will 
scale as vL v = (Mo)vLu. This latter scaling of course is to 
be expected, since the total emitted energy rate can only 
depend on the matter accretion rate, i.e., on the available 
energy. Therefore, the integrated luminosity can be consider 
as a measure of the accretion rate (independent of the mass 
of the central object). 

Fig. 14 shows some indicative spectra of accretion discs 
that are face-on with respect to the observer and extend 
from some inner radius up to some outer radius. The ver¬ 
tical axis is the logarithm of the integrated luminosity in 
units of eV 2 while the horizontal axis is the logarithm of the 
mass-reduced energy of the emitted photons as measured by 
the observer in eV km 1 / 2 . The models that we have used to 
construct the spectra are the same as the models that were 
used to produce the temperature profiles in Fig. 11. The two 
top plots have spin and quadrupolar deformation that cor¬ 
respond to models which have an ISCO, while for the two 
bottom plots the disc terminates on the surface of the star. 
The outer radius of the disc is at around 1000M. In the spec¬ 
tra, three typical regions stand out. At the low energy side 
we can see the low energy tail of the blackbody spectrum 
from the region of the outer radius of the disc, which should 
scale as ~ e 3 (the slope in the plot is ~ 3). At the high en¬ 
ergy side, we can see the high energy exponential behaviour 
of the inner edge of the disc. Finally, between the two we 
can see a modified blackbody spectrum that follows a power 
law ~ s a . The high energy side of the spectrum traces the 
behaviour of the maximum temperature of the disc, which 
is higher the deeper in the gravitational potential the inner 
edge of the disc reaches. Also, the peak of the integrated 
luminosity is consistent with the assumed accretion rate of 
Mo = 10 _12 Mgyear _1 . 

From these spectra we could obtain two observables. 
The energy of the photons where the maximum of the 
integrated luminosity is observed and the value of that 
maximum integrated luminosity. These quantities are given 
in contour plots in terms of the spin parameter and the 
quadrupolar deformability in Figs 15 and 16 respectively. 

Specifically, Fig. 15 gives contours of the mass-reduced 
energy of the photons at the maximum of the spectrum in 
units of eV km 1 / 2 , while Fig. 16 gives contours of the max¬ 
imum integrated luminosity divided by the mass accretion 
rate, uLu/Mo- This last quantity is essentially an alternative 
measure of the efficiency of the accretion, which we saw in 
section 3.6. 

The analysis performed here and the spectra obtained 
demonstrate that one could have an analytic parametrised 
description of an accretion disc around a NS in a way that 
is independent of the EoS and depends only on the NS’s 
spin and quadrupole. Having thus an analytic spacetime that 
describes the geometry exterior to NSs is a very useful tool. 
One could further use this spacetime to do more detailed 
and realistic calculations using ray-tracing and taking into 
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Figure 14. Indicative spectra produced for the same models as the ones used to produce the temperature profiles in Fig. 11. The plots 
show the logarithm of the integrated luminosity i/L„ in eV 2 for a given accretion rate, against the mass-reduced observed photon energy 
Ehv in eV km 1 / 2 . The accretion rate used is Mo = 10“ 12 Algyear —1 and the inclination of the disc with respect to the observer is assumed 
to be zero (face-on). We should also note that leV = 1.6 X 10 12 erg and that leV/h = 2.4 X 10 14 .s 1 , where h is Planck’s constant. 


account all the other effects that we have neglected here for 
simplicity. 

Expressing the various observables of this and the previ¬ 
ous sections in terms of a two dimensional parameter space 
can be quite powerful in probing the properties of NSs in 
GR. We will try to demonstrate this in the following sec¬ 
tion. 


4 COMBINING OBSERVABLES TO 
CONSTRAIN THE EOS 

In section 3 we discussed various possible observables associ¬ 
ated with the orbital motion around a NS, in conjunction to 
some NS properties discussed in the appendix. The observ¬ 
ables/properties are, the location of the ISCO, Risco, the lo¬ 
cation of the surface of the star on the equatorial plane, i.e., 
the stellar equatorial radius R eq , the orbital and the nodal 
precession frequencies of the innermost available circular or¬ 
bit, which depending on the model is either the location of 
the ISCO or the stellar surface, the rotation frequency of 
the star, the orbital and periastron precession frequencies at 
the location where the nodal precession frequency becomes 
zero, the efficiency of a thin accretion disc around the star, 


the maximum effective temperature of the disc, the photon 
energy at the maximum of the integrated luminosity and 
finally the maximum integrated luminosity itself. All these 
quantities were mass-reduced, therefore for our analysis we 
have assumed throughout that the mass of the star is inde¬ 
pendently known. 

But even after having an independent measurement of 
the mass, not all of these quantities are immediately ob¬ 
servable, as is the case for the Risco for example. Some 
of these though can be indirectly related to astrophysical 
observables. The various frequencies for example could be 
observable if they are associated with QPOs, as is assumed 
in the context of the relativistic precession model (Stella 
(2001)). Correspondingly, measurements of the X-ray spec¬ 
tra of LMXBs could provide the efficiency and the integrated 
luminosities or the relevant photon energies. Furthermore it 
might be possible in the future to perform disc tomography 
and obtain thus the temperature distributions of the discs 
in LMXB systems. 

Observing the figures in section 3 one can see that there 
are quantities that we will call here “conjugate”, i.e., they 
give contours that seem to be almost parallel, while there 
are other quantities that we will call “orthogonal”, i.e., their 
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Figure 15. Contours of the mass-reduced photon energy in eV 
km 1 ' 2 where the maximum energy output is observed in the emit¬ 
ted spectrum (which is essentially the photon energy where we 
have the maximum of the spectrum). As before, we assume 
an accretion rate of Mq = 10 —12 Mgyear" 1 . 



Figure 16. Contours of maximum uL^/Mo, expressed in per¬ 
centile, from the spectrum. This quantity is essentially an alter¬ 
native measure of the efficiency of the accretion. 

contours have clear intersections and form a grid. These 
properties of the observables can be used to perform consis¬ 
tency checks and “measurements”. The quantities that are 
conjugate can be used for consistency checks that could dis¬ 
tinguish between proper identifications of the astrophysical 
observables to the corresponding geometric quantities on the 
one hand and possible misidentifications on the other hand. 
For example, the observed QPO frequencies could be caused 
by a different effect than relativistic precession. In that case, 
the comparison between conjugate quantities could help to 
realize the misidentification. In Fig. 17 we present some of 
these conjugate quantities discussed in the previous section. 
We present these plots only for illustrative purposes and for 


that reason as well as for clarity of presentation, we have 
removed the contour values (which can be found in the fig¬ 
ures in the previous section). In the plot on the left of Fig. 
17 we can see the various quantities that are related to the 
radiation emitted from the accretion disc. The plot in the 
middle is the same figure as Fig. 9 and it shows the or¬ 
bital frequency and the periastron precession frequency at 
the location where the nodal precession becomes zero. The 
plot on the right shows the contours of the orbital frequency 
of the innermost possible orbit, against the contours of the 
maximum effective temperature of the accretion disc. 

We should emphasize again that these conjugate prop¬ 
erties could be used as a selection tool. It is important to 
make sure that the quantities that we have measured corre¬ 
spond to the theoretical quantities that are calculated here. 
In a different case, any analysis based on misidentified ob¬ 
servables would give results that would not correspond to 
the properties of the observed system. In the literature of 
the continuum fitting method, there are some selection cri¬ 
teria that have to do with the assumption of a thin accretion 
disc. The main requirement is that the luminosity should be 
L/Le < 0.3, where Le is the Eddington luminosity, so that 
the thin disc assumption does not break down (see the re¬ 
view by McClintock et al. (2014) and references therein). In 
our case, we can use these different conjugate quantities as 
an additional criterion on whether the observed properties 
of a particular system are those that we assume in the the¬ 
oretical analysis presented here. For example, if the QPO 
frequency that we measure from a system is not the orbital 
frequency of the innermost circular orbit, then we would 
not expect it to have a correlation to the maximum effective 
temperature of the accretion disc as the one that we observe 
in the right plot of Fig. 17. 

On the other hand, orthogonal quantities can be used to 
perform measurements of the model parameters. The mea¬ 
surement of two such quantities would give an intersection 
point in the parameter space of (j, a) and therefore a mea¬ 
surement of the j and a of the particular NS. The additional 
knowledge of the mass of the NS would imply the knowledge 
of the first three multipole moments of the star, i.e., the mass 
M, the angular momentum J = jM 2 and the quadrupole 
Q = M 2 = —aj 2 M 3 . The knowledge of these three pa¬ 
rameters for a large enough number of NSs can result in 
constraints for the EoS (Pappas & Apostolatos (2014)). We 
should note that this approach on calculating the first three 
multipole moments of the spacetime and of the central com¬ 
pact object, is an alternative to the approach proposed by 
Pappas (2012) for fitting the evolution of QPO frequencies, 
that was further explored by Boshkayev et al. (2014) and 
Boshkayev et al. (2015). 

In Fig. 18 we show some combinations of these orthogo¬ 
nal quantities that could be used to perform measurements. 
The plot on the left shows contours of the orbital and the 
nodal precession frequencies of the innermost orbit around 
the star and of the rotational frequency of the star. If for a 
particular system we were to identify for example three QPO 
frequencies that would correspond to this triplet, then the 
intersection point of the three corresponding curves would 
indicate a pair of ( j,a ). An alternative combination is pre¬ 
sented in the plot on the right, where we show contours of 
the nodal precession frequency of the innermost orbit around 
the star, contours of the rotational frequency of the star and 
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Figure 17. Left: Contour plots of the accretion efficiency rj (dotted lines), the maximum effective temperature of the disc T^ ax (solid 
lines) and the maximum integrated luminosity vL v /Mq (dashed lines). This plot is a combination of the plots in Figs 10, 12, and 16. 
One can see that these quantities are in a sense conjugate, since their contours are almost parallel. Middle: Plot of the orbital frequency 
(dashed black curve) and the periastron precession frequency (solid red curve) at the location where the nodal precession becomes zero. 
This is the same plot as in Fig. 9. Right: This plot shows the maximum effective temperature of the disc (solid lines) and the orbital 

frequency i/^ ax (dashed lines) of the nearest possible stelar orbit. These two quantities are not as clearly conjugate as the quantities 
in the other plots, but they are clearly not orthogonal. All the quantities plotted and compared are appropriately mass-reduced. The 
numerical values of the relevant quantities on each curve have been omitted for clarity of presentation, but can be found in the figures 
of section 3. 




Figure 18. Left: Contour plots of the orbital frequency at the orbit closest to the stellar surface (dashed black lines), the nodal precession 
frequency at the same orbit (solid black lines), and the rotation frequency of the star itself (dotted red lines). Right: Contour plots of 
the maximum integrated luminosity (solid black lines), the nodal precession frequency at the orbit closest to the stelar surface (dashed 
black lines), and the rotation frequency of the star itself (dotted red lines). These quantities are orthogonal to each other and allow for 
the possibility of inferring j and a. All the quantities plotted and compared are appropriately mass-reduced. The numerical values of 
the relevant quantities on each curve have been omitted for clarity of presentation, but can be found in the figures of section 3. In each 
plot we represent with three thick curves the curves that correspond to some measurement of each of the quantities. The intersection of 
these three curves indicates a measurement of j and a. 


contours of the maximum integrated luminosity. Again we 
can see that a possible measurement of these three quantities 
could suggest a pair of ( j , a). 

The choice of the quantities plotted in Fig. 18 is only 
one possible combination. For example, one could choose 
to plot the orbital or the precession frequencies at different 
locations, depending on the preferred model for the QPO 
frequencies (for an alternative, see for example the work by 
Stuchlik et al. (2015)). This demonstrates the power of hav¬ 


ing a parameterized analytic description of the spacetime, 
since one could have as an extra free parameter the location 
at which a QPO is created (assuming the relativistic preces¬ 
sion model for example) and try to fit for that parameter as 
well. 

Another alternative application could be the following, 
if one were to have two well identified quantities, such as 
the stellar rotation frequency and the orbital frequency at 
the innermost circular orbit for example, then one could at- 
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tempt to test a hypothesis about the nature of a third fre¬ 
quency by seeing if the right curve passes through the point 
were the two other curves intersect. This would be an even 
stronger test of the assumptions that enter the modelling of 
the various astrophysical observables, than the tests one can 
do using the conjugate quantities. 


5 CONCLUSIONS 

In the current work we have performed an investigation of 
the astrophysical properties of NSs that are related to the 
spacetime geometry around them. These properties were 
presented separated in two categories, the first being the 
properties of the geodesic motion of a test particle or a fluid 
element around the NS and the second being the proper¬ 
ties of matter accreting onto the NS while emitting thermal 
radiation. 

To calculate the properties of the geodesic orbits and 
of the accretion disc we have used the analytic two-soliton 
spacetime, which was shown by Pappas & Apostolatos 
(2013) to accurately reproduce the spacetime around NSs if 
one chooses the parameters of the two-soliton spacetime in 
such a way so that the spacetime has the same first four non¬ 
zero multipole moments as the corresponding NS. Further¬ 
more, instead of assigning specific values for the multipole 
moments of the NS, we have taken advantage of the recently 
found universal relations between the multipole moments 
(see work by Pappas & Apostolatos (2014) and Yagi et al. 
(2014)). By further mass-reducing all the observable quan¬ 
tities we have managed to give a description of the various 
orbital and accretion properties in an EoS-independent way 
that depends on only two parameters, the spin parameter of 
the NS, j. and the quadrupolar deformability, a. 

Since the various observables are parameterized only by 
the spin parameter and the quadrupolar deformability, one 
can construct contour plots for these observables on the pa¬ 
rameter space of (j. a). The utility of these contour plots is 
twofold. The first is that one could use these plots to test 
the assumptions that enter the modelling of the various as¬ 
trophysical observables, such as that QPOs originate from 
geodesic motion for example or assumptions about the prop¬ 
erties of accretion discs. The second is that by combining the 
right observables, one could in principle measure the spin 
parameter and the quadrupolar deformability (from the in¬ 
tersection of the curves of two or more observables that have 
been measured). This application assumes the independent 
knowledge of the mass of the NS. Therefore once there ex¬ 
ists an independent mass measurement, the measurement of 
j and a implies the knowledge of the first three multipole 
moments of the NS, i.e., the mass M, the angular momen¬ 
tum J = jM 2 , and the quadrupole moment Q = —aj 2 M 3 . 
The measurement of the first three moments of a variety of 
NSs will be able to constrain the EoS for the matter inside 
NSs. 

Furthermore, measuring j and a will give us information 
for other properties of NSs that are not directly measurable, 
such as the radius of the NS (see appendix A), which enter 
the calculations related to other astrophysical processes such 
as photospheric radius expansions. This would provide the 
possibility of cross-checking and refining our models. 

All these possibilities, and possibly even more than 


those this author can conceive, are available to us because 
the analytic two-soliton spacetime and the universal rela¬ 
tions between the multipole moments allow for an analytic 
parametric description of the exterior of NSs in an EoS- 
independent way with a minimal set of parameters. The re¬ 
sults presented here could be extended further by extending 
the analysis and the modelling of the accretion disc to in¬ 
clude more phenomenology, such as for example taking into 
account temperature colour corrections or ray tracing for the 
photons and so on and so forth. 

Finally, this analysis could probably be extended to 
quark stars as well. On the one hand, it is likely that the 
two-soliton spacetime could describe the exterior of quark 
stars, if the right value for the multipole moments is cho¬ 
sen (although this would need to be demonstrated). On the 
other hand, the universal relations for the multipole mo¬ 
ments, used here for NSs, were shown by Yagi et al. (2014) 
to hold for quark stars as well (with a slight modification). 
Therefore, one could perform a similar analysis to include 
quark stars. The difference with respect to the analysis pre¬ 
sented here would be on the behaviour of the surface of quark 
stars, which one would expect it to be at smaller radii than 
it is for NSs. This means that maybe the surface of a quark 
star would not interfere with the corresponding ISCO, which 
would simplify the analysis. 
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APPENDIX A: EQUATORIAL RADII OF NSS 

The presentation in this work is based on the fact that it 
is possible to have an EoS-independent description of the 
spacetime around a NS. Some of the results presented here 
also depend on whether the radius of the star is such that 
it does not hide the region of interest. Therefore it would 
be useful if we could have an EoS-independent parameteri¬ 
zation of the radius of a NS (in particular the reduced cir¬ 
cumferential radius R c irc/M) in terms of the spin parameter 



Figure Al. The plot on the left shows contours of the fit of the 
normalized circumferential equatorial radius R c i rc /M of NSs as a 
function of the spin parameter j and the quadrupolar deformation 
parameter a. The plot on the right shows overlaid to the normal¬ 
ized radius (dotted contours), contours of relative difference per 
cent between the actual radius and the fit (solid contours). One 
can see that the accuracy of the fit is almost for the entire pa¬ 
rameter space better than 6 per cent with a vertical region on the 
side close to maximal rotation that is between 6 and 10 per cent. 


j and the quadrupolar deformability parameter a. It turns 
out that such a parameterization is possible. 

From the set of NS models that were con¬ 
structed numerically using realistic EoSs in the work by 
Pappas & Apostolatos (2014) (using the RNS numerical 
code developed by Stergioulas & Friedman (1995)), we plot¬ 
ted the reduced radii as functions of the spin parameter j 
and the square root of the quadrupolar deformability pa¬ 
rameter y/a. All the models from all the EoSs appear to 
approximately occupy the same simple surface. Therefore 
we have fitted the surface with a function of the form, 

R sq /M = Bo + Bij + B 2 j 2 + (A + Aij + A2?){.y/af ri 

+(Co + Cij + C 2 /)(^)V (Al) 

where the values of the parameters are, Ao = 0.00927584, 
Ai = -0.0252801, A 2 = 0.0497335, B 0 = -0.358824, Bi = 
3.15892, B 2 = -5.30171, C 0 = 2.94923, Ci = -3.20369, 
C 2 = 6.02522, Afi = 4.12566, and JV 2 = 0.996284. For these 
values, the accuracy with which the surface reproduces the 
radii of the NSs is better than 10 per cent, with most of 
the parameter space being described with accuracy better 
than 6 per cent. It is possible that with a better choice of a 
fitting function better accuracy can be achieved, but for the 
purposes of this work the current accuracy will be adequate. 
The resulting distribution of NS equatorial radii with respect 
to j and a can be seen in the plot on the left of Fig. Al. On 
the right we have for illustrative purposes a contour plot of 
the relative difference between the fit and the actual radii 
values. The contours indicate the regions where the error is 
between 6 and 10 per cent. 

This parameterization of the equatorial radius could be 
also seen in the spirit of reducing the number of free pa¬ 
rameters in describing NSs, as it was the aim of the work 
by Baubock et al. (2013). There has been a similar attempt 
in expressing the NS radii by Chakrabarti et al. (2014) in 
an EoS-independent way, the difference being that the ra¬ 
dius was given as a function of the spin parameter and the 
rotation frequency of the star. 
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Figure Bl. Same figure as Fig. A1 but for the stellar rotation 
frequency. The plot on the left shows contours of the normalized 
rotation frequency M X /(km X kHz) of NSs as a function of the 
spin parameter j and the quadrupolar deformation parameter a. 
The plot on the right shows overlaid to the normalized rotation 
frequency (dotted contours), contours of relative difference per 
cent between the actual rotation and the fit (solid contours). One 
can see that the accuracy of the fit is for the entire parameter 
space better than 1 per cent. The plot shows only contours of 
relative difference larger than 0.7 per cent. 




APPENDIX B: ROTATION FREQUENCY OF 
NSS 

As it was necessary to have an EoS-independent description 
of the circumferential radius of the star, it is also neces¬ 
sary to have an EoS-independent parameterization of the 
stellar rotation rate in terms of the spin parameter j and 
the quadrupolar deformability parameter a. As for the case 
of the surface radius, from the set of NS models that were 
constructed numerically using realistic EoSs in the work 
by Pappas & Apostolatos (2014), we plotted the normalized 
frequency M x f/j, where M is the mass of the star, / is the 
rotation frequency, and j is the spin parameter, as functions 
of the spin parameter j and the square root of the quadrupo¬ 
lar deformability parameter y/a. All the models from all the 
EoSs appear to approximately occupy the same simple sur¬ 
face and we have fitted the surface with the same function 
as the one we used for the radius, which is of the form, 

= *8 0 + VSij + © 2 / 2 + (2l 0 + 2lij + 2l2j 2 )(%/a) ni 

J 

+(£o + tij + £ 2 j 2 )(Va) n2 . (Bl) 

The values of the parameters of the fit are, 2lo = 15.0297, 
2li = -0.114154, 2h = -7.72439, ©o = -1.48338, «i = 
-1.07874, <8 2 = 1.64592, £ 0 = -2.45303, £i = 3.40995, 
£2 = 7.39354, m = —1.1698, and 112 = —4.50216. For these 
values for the fit parameters, the accuracy with which the 
surface reproduces the normalized frequency M x f/j of the 
NSs is everywhere better than 2 per cent and in particular, 
in the region of the parameter space that is of interest to 
our analysis, it is better than 1 per cent. The resulting dis¬ 
tribution of NS equatorial radii with respect to j and a can 
be seen in Fig. Bl. 

One could note here that oc 1/1, where I is 

the reduced moment of inertia as defined by Yagi & Yunes 
(2013b, a). 


APPENDIX C: NODAL PRECESSION 
FREQUENCY Q. z AROUND NSS 

It was found by Pappas (2012) (and has been recently redis¬ 
covered by Gondek-Rosinska et al. (2014); Wisniewicz et al. 
(2015)) that the significant deviation of the quadrupole of a 
NS from the quadrupole of a Kerr black hole with the same 
mass and angular momentum, can result in a very differ¬ 
ent behaviour of the nodal precession frequency from what 
one would observe for a Kerr spacetime. In particular, if the 
quadrupolar deformation of the compact object due to rota¬ 
tion is high enough, this can result in a vertical oscillation 
frequency k z that is higher than the orbital frequency B in 
the region close to the ISCO. Then as we move further away 
from the ISCO, at some point the vertical frequency and the 
orbital frequency become equal and from that point outward 
the orbital frequency becomes higher than the vertical fre¬ 
quency. That means that the nodal precession, = O — k z , 
in the outer region of the spacetime is towards one direction, 
at the point where the two frequencies are equal there is no 
precession, and in the inner region the precession is in the 
opposite direction. 

This effect can be of great significance in modelling ac¬ 
cretion discs and would constitute a very clear signature of 
geodesic precession if it were to be observed in QPOs, as 
it was concluded by Pappas (2012). For the sake of com¬ 
pleteness and in order to make the effect more clear so as 
to facilitate better understanding of the results presented in 
the main text, we will give here a brief description of the be¬ 
haviour of the geodesic precession frequencies for different 
rotation rates and different quadrupolar deformations. 

Figs Cl and C2 present the orbital frequency and the 
precession frequencies as functions of the distance from the 
central object (top panels of each figure) as well as the two 
precession frequencies as functions of the orbital frequency 
(bottom panels of each figure). All frequencies and distances 
have been normalized with respect to the mass M of the 
compact object which provides a natural length scale in ge¬ 
ometric units (assumed to be in km). The different frequen¬ 
cies shown here are related to the angular frequencies as 
B = 2ni/. The item of interest in all figures are the curves 
for the nodal precession frequency, v z = v v — k z /(2n), while 
the other frequencies are shown for illustrative purposes and 
perspective. 

Fig. Cl shows how the precession frequencies behave for 
two different spin parameters ( j = 0.3 left column, j = 0.4 
right column) if we change the quadrupolar deformation pa¬ 
rameter a. The values for a are in the range [1.2,10], where 
the value a = 1 would correspond to the Kerr quadrupole 
M 2 ~ —j 2 M 3 . We have chosen these values for the spin 
parameter because they don’t correspond to extremely 
rapidly rotating NSs. In particular, j = 0.4 is well within 
the region where the zero of the nodal precession is outside 
the surface of the star, as Fig. 9 shows. We should also note 
here that lower values of a correspond to NSs with masses 
that are towards the maximum mass/high density end of 
the mass-radius diagram, while higher values of a corre¬ 
spond to lower mass NSs (see Laarakkers & Poisson (1999); 
Pappas & Apostolatos ( 2012 ); Pappas & Apostolatos 
(2012); Urbanec et al. (2013)). 

Returning to the behaviour of the nodal precession, we 
can see from Fig. Cl that for a given spin parameter, for 
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Figure Cl. The figure shows the various geodesic normalized frequencies for NS spacetimes. The two top panels show the orbital 
frequency (red), the periastron precession frequency Mv p (blue), and the nodal precession frequency Mv z (green) as functions of 

the normalized circumferential distance R c i rc /M from the central object. In the case of the nodal precession, since in some regions it 
becomes negative, we have plotted the absolute value. The two bottom panels show the two precession frequencies (same corresponding 
colours) as functions of the orbital frequency. The spikes in the plots of the nodal precession indicate where the precession becomes zero. 
The figures in the left column correspond to compact objects with the same value of the spin parameter j = 0.3, while the figures on 
the right column have j = 0.4. The different curves in each figure correspond to different values of the quadrupolar deformability in the 
range 1.2 ^ a. ^ 10. 





Mxv^p (kmxHz) 


Mxv^p (kmxHz) 


Figure C2. Same figure as figure Cl, with the difference being that the figures in the left column correspond to compact objects with 
the same value of the quadrupolar deformability a = 5, while the figures on the right column have a = 9. The different curves in each 
figure correspond to different values of the spin parameter in the range 0.05 ^ j ^ 0.75. 
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small values of the deformation parameter, the precession 
frequency behaves as a power law just like the nodal pre¬ 
cession frequency of the Kerr geometry (the leading order 
behaviour is of the form ~ r~ 3 for Kerr). As a increases 
though, this power law starts to break at the ISCO (which 
is at the smallest radius side of the curves, or highest v v ) 
until for some value of a the first spike appears at a radius 
outside the ISCO, which corresponds to a zero nodal pre¬ 
cession. Beyond that point, for even higher values of a the 
position of the zero of the nodal precession moves further 
out until we reach the maximum value of a. One should also 
note that at the far away region (larger radii, smaller or¬ 
bital frequencies) the curves seem to bundle together, which 
is what one would expect since at that region the frequencies 
are mainly characterized by the rotation and the quadrupole 
or any higher moment should be less important. 

Fig. C2 shows again the various frequencies but in this 
case we have plotted for two values of the quadrupolar de- 
formability parameter, a = 5 (left column) and a = 9 (right 
column), models with varying spin parameter in the range 
j £ [0.05, 0.75]. To give some idea of what sort of NSs that 
would correspond to, depending on the EoS, these values of 
a correspond to NSs with masses between 1 — 1.4Mq. 

One can notice (particularly from the far away be¬ 
haviour of the nodal precession frequencies) that the first 
effect of increasing rotation is to move upwards in a paral¬ 
lel fashion the frequency curves. This is something that was 
expected from the behaviour of the Kerr nodal precession 
frequencies, since the spin parameter multiplies the lead¬ 
ing order behaviour of the frequencies (~ jr~ 3 ). As in the 
previous case, we notice that as the rotation increases and 
with it the magnitude of the quadrupole moment (which is 
oc j 2 ), the power law behaviour develops a break close to 
the ISCO, which then propagates outwards until a zero of 
the precession frequency appears (in the form of a spike in 
these log-log plots) which moves outwards, to larger radii 
or lower orbital frequencies, as the spin parameter continues 
to increase. The difference between the two columns with 
the different values of a is that in the case of the higher 
quadrupolar deformability (right column plots) the position 
where the precession becomes zero can be further out than 
in the case of the lower a. 

Finally for actual NSs, there is also the issue of the 
location of the surface of the star (presented in Appendix 
A). It can be for some EoSs that the NS that corresponds 
to a combination of a spin parameter and a quadrupolar 
deformation that lead to a vanishing precession frequency at 
some radius, has actually an equatorial radius that is further 
out. In these cases, the geodesics, the precession frequencies 
of which we are investigating, terminate at the surface of the 
star. 

This paper has been typeset from a IpX/r^TpX file prepared by 
the author. 
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